Induced Triplet Pairing in clean s-wave Superconductor/Ferromagnet layered 

structures 
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We study induced triplet pairing correlations in clean ferromagnet/superconductor/ferromagnet 
heterostructures. The pairing state in the superconductor is the conventional singlet s-wave, and 
the angle a between the magnetizations of the two ferromagnetic layers is arbitrary. We use a 
numerical fully self-consistent solution of the microscopic equations and obtain the time-dependent 
triplet correlations via the Heisenberg equations of motion. We find that in addition to the usual 
singlet correlations, triplet correlations, odd in time as required by the Pauli principle, are induced 
in both the ferromagnets and the superconductor. These time-dependent correlations are largest 
at times of order of the inverse of the Debye cutoff frequency, u>d, and we find that within that 
time scale they are often spatially very long ranged. We discuss the behavior of the characteristic 
penetration lengths that describe these triplet correlations. We also find that the ferromagnets 
can locally magnetize the superconductor near the interface, and that the local magnetization then 
undergoes strongly damped oscillations. The local density of states exhibits a variety of energy 
signatures, which we discuss, as a function of ferromagnetic strength and a. 

PACS numbers: 74.45.+C, 74.25.Bt, 74.78.Fk 



I. INTRODUCTION 

Triplet Cooper pairing is no new phenomenon: it has 
long been recognized to be responsible for superfluidity 1 
in 3 He as well as for superconductivity in some electronic 
materials. This occurs when the pairing interaction is in 
a partial wave with odd £. However, recent observations 
have raised the possibility of induced triplet pairing cor- 
relations in s-wave superconductors. It is a matter of 
elementary physics that the Cooper pair wavefunction 
must be antisymmetric under exchange of the two elec- 
trons to satisfy the Pauli principle. For spatially symmet- 
ric s-wave superconductors with decoupled spatial and 
spin degrees of freedom, the spin singlet pair is the only 
possible antisymmetric state of an electron pair. Triplet 
pairing states on the other hand, where the spin state 
is symmetric, are obviously allowed when the pairing is 
spatially antisymmetric, such as in p-wave superconduc- 
tors. Triplet states in systems with s-wave pairing, even 
in momentum or coordinate space, would naively appear 
to violate the Pauli principle. However, many years ago 
Berezinskii proposed 2 - a triplet state in superfluid 3 Hc, 
which involved spatially symmetric correlations. Berezin- 
skii's triplet pairing correlations, involving different-time 
pairing, did not violate the Pauli principle by virtue of 
being odd in time, thus allowing a triplet state in a sys- 
tem with s-wave interactions. While such a state did not 
turn out to be appropriate to describe superfluidity in 
3 He, its consideration has led the way to the study of 
cases where some sort of time-reversal symmetry break- 
ing mechanism may allow an odd time triplet state to be 
induced in systems with spatially symmetric interactions. 

Interest in exotic triplet pairing arises from many quar- 
ters. In the case of two component cold atomic gases 3 



with short-range s-wave interactions, in which the two 
species have the same mass but different chemical po- 
tentials, it may be possible to induce a triplet pairing 
that breaks time reversal symmetry. In electronic materi- 
als, an important issue is the possible existence of a long 
range proximity effect in Superconductor/Ferromagnet 
(SF) heterostructures. Interest in these heterostructures 
arises in turn from their possible* applications. Of par- 
ticular interest is that the thermodynamic and transport 
properties of FSF trilayers are found to depend strongly 
on the relative orientation of the magnetization in the 
two F layersi^i^ This rather well-understood ^ 10 ' 11 ! 12 
fact makes these structures candidates as spin valves. 

There have been no unambiguous observations of in- 
duced triplet correlations in SF structures involving s- 
wave superconductors. However, there have been some 
enticing experimental hints in the form of long ranged 
proximity induced superconducting behavior in SF mul- 
tilayers with strong exchange fields. The observed ef- 
fects are over length scales much larger than those of 
the usual SF proximity effect and more like the much 
longer length scales associated with the standard prox- 
imity effect between a superconductor and a normal non- 
magnetic metal. These observations include measure- 
ments in superlatticesi 3 - with ferromagnetic spacers and 
SQUIDs 14 with ferromagnetic interlayers. Superconduct- 
ing characteristics, such as a critical temperature and en- 
hanced sub-gap conductance, have been observed in point 
contact conductance measurements on s-wave supercon- 
ductor/half metallic systems Perhaps most compelling 
is the observation of a Josephson current through a strong 
almost half-metallic, ferromagnet.— All of these experi- 
ments indicate long range superconducting correlations 
that are not destroyed by a strong exchange field: a 



2 



triplet state would obviously be consistent with these ex- 
periments. To fully understand the behavior of these 
systems, further studies involving, for example, quanti- 
ties sensitive to the gap such as the local density of states 
(local DOS), are needed. Until both better theoretical 
models and more varied experimental observations are 
made available, one cannot conclusively say that there is 
indeed an induced triplet state in these systems, but at 
this point the facts fit this explanation and no better one 
has been proposed^ 

Many theoretical studies agree that it is pos- 
sible to induce this exotic state in certain SF 
systemsii^^^^^^^^^S^S (perhaps even in 
the nonmagnetic^ case) with ordinary singlet pairing in 
S. Some studies use an SFF'S arrangement in which the 
F and F ' layers have different magnetization orientations. 
Others assume that a domain structure in a single F layer 
is responsible for the symmetry breaking. Yet others as- 
sume an FSF system with different in-plane magnetiza- 
tion orientations in the F layers. Whatever the mecha- 
nism for the symmetry breaking, such arrangements can 
induce, via proximity effects, triplet correlations of dif- 
ferent kinds. Recently it was shown that a Josephson 
supercurrent can exist in a half metal by virtue of equal 
spin triplet pairs and spin flip scattering events at the 
interfaces^ To understand and probe the underlying 
triplet state, investigations have been done on conduc- 
tance spectra in simpler FS structures^ 1 with arbitrary 
magnetization alignment and with spin active interfaces, 
as well as in diffusive^ SF junctions, through character- 
izing the possible superconductor symmetry classes con- 
sistent with Pauli's principle. With the exception of our 
early worki£ on SFS trilayers and some recent work on 
SFS Josephson junctions^ the studies above are done 
in the dirty limit through linearized Usadel-type or other 
quasiclassical equations. The disadvantage of a quasiclas- 
sical approach is that it is unsuitable for magnets with an 
exchange field on the order of the Fermi energy. Thus, it 
cannot properly model a strong ferromagnet, and it does 
not allow for atomic scale oscillations in the pair ampli- 
tude. A good quantitative explanation requires a self- 
consistent treatment of a fully microscopic model. Thus, 
as pointed out in a recent review^ the very existence of 
triplet correlations in clean FS structures was until very 
recently generally doubted, and these doubts have only 
very recently^ been dispelled. 

In this paper we explore the phenomenon of induced 
triplet correlations, odd in time, of clean FSF structures, 
where S is an ordinary s-wave superconductor and the 
magnetizations in the two F layers are rotated by an arbi- 
trary angle a. We assume strong ferromagnets (up to the 
half-metallic limit), and smooth, sharp interfaces. In this 
geometry, triplet correlations with total spin projection 
m = on the axis of quantization of the Cooper pairs are 
in general possible and, when the relative magnetizations 
(which we assume as usual are both parallel to the inter- 
faces) are not aligned, triplet components with m = ±1 



are allowed also. To satisfy the Pauli principle, these 
spatially symmetric triplet pairing correlations must be 
odd in frequency or time 3, That such correlations are al- 
lowed, does not mean that they must exist, nor that they 
must exist over an extended spatial range. We find, how- 
ever, via a fully self-consistent solution to the microscopic 
Bogoliubov de-Gennes^ (BdG) equations that such cor- 
relations do indeed exist, and that the penetration depth 
associated with them can be very long. Our use of the 
BdG equations allows us to study strong ferromagnets. 
Self-consistency is fundamental: non-self-consistent solu- 
tions are found to violate the Pauli principle. Thus, the 
time consuming step of calculating fully self-consistent 
solutions is necessary to properly model the proximity ef- 
fects which allow for the mixing of superconducting and 
ferromagnetic orderings that causes these induced corre- 
lations. 

In Sec. [IT] of this paper, we discuss the basic equations 
and our method for numerical self-consistent solution. 
There, the extraction of the all-important time depen- 
dence via solution of the Heisenberg equations of motion 
for the relevant operators is explained in detail. Expres- 
sions for all of the time-dependent triplet correlations are 
also derived. The equations for the local density of states 
(local DOS) and the local magnetic moment (which we 
use to discuss the reverse proximity effect, that is, the 
penetration of the magnetism into the superconductor) 
are also presented. In the next section (Sec. IIII[) we be- 
gin by presenting an extensive discussion of the triplet 
correlations as a function of position, a, and magnet 
strength. The appropriate penetration depths are ex- 
tracted and discussed. Results for the DOS, the magnetic 
moment and the temperature dependence of both triplet 
and ordinary singlet correlations are also given. Finally, 
in Sec. IIV[ a brief conclusion and summary is given. 



II. METHODS 

The geometry we consider consists of a planar FSF 
junction as depicted in Fig. [TJ The thickness of the su- 
perconducting layer is ds and the F layers have thick- 
nesses dp 1 and dj? 2 . The system is assumed to be infinite 
in the plane perpendicular to the layers, which we label 
as our x — z plane. The magnetizations of the F layers, 
which are in this plane, form an angle ±a/2 with the z 
axis, which is that of the direction of quantization of the 
spins. 

Our starting point is the Bogoliubov de-Gennes (BdG) 
equations^ for the system under consideration. The 
derivation of the BdG equations for the case of interest 
requires some care with the conventions for all operator 
phase factors, which are not universally agreed upon in 
the literature, and which may give rise to different signs 
in some of the equations below. We write the effective 
BCS Hamiltonian, 7Y e ff, as 
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FIG. 1: (Color online) Schematic of the FSF junction. The 
y axis is normal to the interfaces. The left ferromagnet layer 
denoted Fi has a magnetization oriented at an angle —a/2 in 
the x — z plane, while the other magnet F2, has a magnetiza- 
tion orientation at an angle a/2 in the x — z plane. All layer 
widths are labeled. 



^eff = / <?r{£^(r)Wa(r) + iEK^^ (2.1) 

a a, (3 a,/3 
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where H c = -l/(2m)V 2 - E F + U(r), a are the set 
of Pauli matrices, spin is denoted by Greek indices, and 
as usual, we represent the magnetism of the F layers by 
an effective exchange Stoner energy h(r) which will in 
general have components in both the transverse [x, z) 
directions. The spin independent scattering potential is 
denoted U(r), and A(r) is the usual pair potential. 

To diagonalize the effective Hamiltonian, the field op- 
erators tj}' a and ip a are expanded by means of a Bogoli- 
ubov transformation, which, for our phase convention, 
we write^ as: 

( r ) = ( U »T ( r )^« ~ U "T ( v )ln) , ( 2 - 2a ) 

n 

</U M = ( r )^" + v ^ , (2.2b) 



where u na and v na are the quasiparticle and quasiholc 
amplitudes, and 7„ and 7^ are the Bogoliubov quasipar- 
ticle annihilation and creation operators, respectively. 

We require that the transformations in Eqs. (|2.2[) di- 
agonalize 7i e ff, 

[HeS,Jn] = -Cnln, (2.3a) 

[U^ll] = wl (2.3b) 

One can also take the commutator [tp a (r),Ti. c s]- With 
the magnetizations in the x — z plane as explained above, 
this gives the following, 



[VT(r),«eff] = (Wc - M<Mr) - Mi(') + A(r)^{(r), 
[Vi(r),Weff] = (H c + h z )^{v) - h x ^(v) - A(r)^{(r). 
Inserting (|2.2[) into (|2.4[) and using Eqs. (|2.3p yields the general spin-dependent BdG equations, 



/H -h z (y) ~h x {y) 

-h x (y) H + h z (y) A(y) 

A(y) -(n -h z (y)) 

V A(y) -My) 




"«t(2/) n 

"»i(y) 

Wr»T(j/) 



/ "«T(y) N 
"»i(y) 

^n|(2/)/ 



(2.4a) 
(2.4b) 



(2.5) 



where the single particle Hamiltonian Tio is defined as, 

£%). (2.6) 



7Yo = 77^ h £_L — -^F 



2m 



A plane wave factor e 



has been canceled in both 
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sides of Eq. (|2.5p . The longitudinal momentum opera- 
tor, p y , is given by, p y — —id/dy, s± is the kinetic en- 
ergy of the transverse modes, A(y) is the self-consistent 
pair potential, and U(y) is a scalar potential represent- 
ing interface scattering characterized by a delta func- 
tion of strength Hb- The ferromagnetic exchange field 
h(y) = (h x (y),0 7 h z (y)), vanishes in the S layers. We 
have h x (y) = ho sirx(— a/2) and h z (y) — ho cos(— a/2) 
in the Fi layer, where ho is the magnitude of the ex- 
change field, while in F 2 , h x (y) — h s'm(a/2), and 
h z (y) — /i cos(a/2). We refer to Fig. [T] for details. The 
dimensionless parameter I = ho/Ep conveniently char- 
acterizes the strength of the magnetism. One thus has 
/ = 1 in the half metallic limit. If we take a = or ir, 
i.e. the magnetizations of both layers lie along the same 
direction, or if there is only one F layer, then h x = 
and we recover the simpler form of the BdG equations 
used^i in other contexts. We can find the quasiparticle 
amplitudes on a different quantization axis in the x — z 
plane forming an angle a' with z, by performing a spin 
rotation $ n — » U(a')$ n with, 



U(a') = cos(a'/2)i ® 1 - ism(a'/2)p z 



(2.7) 



where we have introduced p as a set of Pauli-like matrices 
in particle-hole space. It is convenient to use the matrices 



p in conjunction with the ordinary Pauli matrices <x in 
spin space to rewrite Eqs. (|2.5j) in the more compact but 
perhaps less transparent way: 

[pz ® (H i - h z a z ) + (A(y)p x - h x l) ® a x ] $„ = e n $„, 

(2.8) 

where $„ = (w„T(y),Wni(y),-iM(y)^nl(y)) T > witn tlie 
superindex denoting transposition. 

The usual self consistency condition relates the spec- 
trum obtained from Eq. (|2.5j) to the inhomogeneous pair 
potential A(y) by an appropriate sum over states: 



A(y) 



ul{y)vi(y)] tanh(e n /2T), 
(2.9) 

where the prime on the sum indicates that only those 
positive energy states with energy less than the pairing 
interaction energy cutoff, u>d, are included, and T is the 
temperature. The function g(y) vanishes in the F layers 
while in the S layers it takes the value of the usual BCS 
singlet coupling constant in the S material. 

With an appropriate choice of basis , 34 > 35 Eqs. (|2.5p can 
be cast into a finite AN x AN dimensional matrix eigen- 
value system. In dimensionless form, it reads, 



J 



Ho — H z —H x 
—H x Ho + H z 
D 
D 




D 

-(H -H z ) 



D 


-H x 
-{H + H z ). 



(2.10) 



where e„ = e„/ Ep, and ^! nj the transpose of 



' U nN' ^ 



nl ' ' 



1 T 
' ' M n7V> v nl> ■ 



T 1 
' ' V nN' V nl> ■ 



J nN> 



(2.11) 



contains the expansion coefficients associated with the set of orthonormal basis functions. We write u"(z) = 
2 / d^2 ± u° q sin(qTT z / d) , and v"(z) = \/2/d^2 „ =1 v" q sin(qnz/d), for a =|, |. The necessary matrix elements 



analogous to Eqn. (|2.5p for different 7r junction geometries and for strictly collinear magnetization orientations have 
been calculated in previous work. 34 ! 35 The situation is more complicated for the case of a FSF trilayer considered 
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here, with the magnetization angles of the two F layers forming an angle a. The matrix elements are then written as, 



(-ffo)mn 



(H z )r 



rmr 
kpd 



5 mn + ZB\Um-n{dF\) + Um-n(dFl + ds) — Um+n{dFl) 



{H x ) r 



(D)r 



pa I Ef 
Um+n(dFl + ds)], 

■■ cos(a/2) \K m - n {dF\) - )C m+n (d F i) + JC m+n (d F i + ds) 
fc m -n(dFi+ds)], m^n, 

dFl + dF2 



t=t" cos(a/2) 
h . 



+ A^2m(d_Fl + ds) — K-2m(dFl) 



m = n, 



rn — n (dFl) + f£ m +n(dFl + ds) 



E F 

K m - n (d F i + d s )], m^n, 

. / , n s \dp2 - dFl 
■ sin(a/2) 



+ K-2m{dFl + d s ) + IC 2m {dFl) 



, m = n, 



E F d 



d F1 +ds 



dy sin 



mixy 



A(y) sin 



nny 



(2.12a) 

(2.12b) 
(2.12c) 

(2.12d) 
(2.12e) 

(2.12f) 



where Zb — 2Hb / '(kpd) is a convenient dimensionless 
measure of interfacial scattering. We have also defined: 



These equations of motion, given Eqns. 



have the 



IC n {y) 



sin 



^)-cos(^) 



(2.13) 



We now consider the appropriate quantities that char- 
acterize the induced triplet correlations. To do this, we 
define the following triplet pair amplitude functions in 
terms of the field operators, 

foM) =i[(^(r,t)^(r,0)) + (^(r,i)^ T (r,0))], 

(2.14a) 

h(r,t) =i[(^(r,t)^(r,0)> - (^M)^(r, 0)>]. 

(2.14b) 

We will later demonstrate that these amplitudes vanish 
at t — 0, as required by the Pauli principle. 

To make use of these expressions, it is most convenient 
to use the Heisenberg picture. Thus we write if) q in the 
Heisenberg representation: 



^(t) = e (^efft)^ e (-iW 6ff t)_ 



(2.15) 



To put this in terms of the quasiparticle amplitudes, we 
apply Eqns. (|2.2[) and the transformation Eqns. (|2.3p . We 
can then immediately write down the Heisenberg equa- 
tions of motion for the 7's as 



'~dT 



off 



and 



(2.16) 



(2.17) 



solutions 7n(i) = r y n e~ zerit and i n (t) = i n e %tnt . When 
we substitute these results into the above equations for 
fo and fi, taking into account Eqns. (|2.2p we obtain: 



f o(y^) =^^2[Unl(y)Vnl(y) - U nl (y)v n i(y)]Cn(t), 
n 

(2.18a) 

h{y,t) = - i^[ u ™T(y) w «T(y) + uni(y)vni(y)](n(t), 

n 

(2.18b) 

where C„(t) = cos(e„t) — i sm(e n t) tanh(e„/2T). The spa- 
tial dependence of the complex quantities to(y,t) and 
fi(y,<) is, in our geometry, on the y coordinate only. 
They vanish identically at t = 0. 

We will focus on in our study of the induced triplet 
correlations on the time dependent quantities fo(y, t) and 
fi(y,t). Their existence at t > is allowed by the Pauli 
principle. It is also important to sort out when it is 
allowed by the spin symmetries: when the axis of quanti- 
zation of the Cooper pairs is the only axis of quantization 
in the system (i.e., when a = 0) then it is not hard to see 
that the total spin operator S of the Cooper pairs does 
not commute with the Hamiltonian. This is best seen 
directly from the matrix expression on the left side of 
Eqn. (|2.5p . On the other hand, S z and the Hamiltonian 
do commute in this case. However, when a (and therefore 
h x ) is nonzero, then no component of S commutes with 
the effective Hamiltonian. From this spin symmetry ar- 
gument it follows that the induced amplitude fi(y, t) may 
exist (at finite times) only at nonzero a, while fo(j/,t) is 
allowed ant any a. For a = tt, when the magnetizations 
are antiparallel and along the x axis, no triplet ampli- 
tudes with nonzero component along that axis can exist. 
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The matrix U (a) in Eq. (|2.7|l can be used to verify this 
by performing the corresponding spin rotations. That 
the existence of certain quantities is consistent with all 
symmetry properties does not mean that these quanti- 
ties will indeed be nonvanishing, and it certainly tells us 
nothing about the possible range and behavior in space 
and time of these amplitudes. To determine this requires 
detailed calculations. 

Also of considerable interest in F/S structures is the 
reverse proximity effect: the leakage of magnetism out 



of the magnets and into the superconductor. This can 
be characterized by the local magnetization m(y). It is 
defined as, 

m = -MflE^H). (2-19) 

where fis is the Bohr magneton. The vector m has two 
components in the FSF geometry discussed. Both com- 
ponents depend on y. They are: 



m z {y) = -vbJ2{ [\ u ni(y)\ 2 - Ki(y)\ 2 ] f n + [k T (y)l 2 - ki(y)l 2 ] (i -/«)}, (2.20) 



and 

m x (y) = -2/i B ^ {i U nl(y)Unl(y)} fn + [Vn^(y)Vnl(y)} (1 - fn)} ■ (2-21) 



It is convenient to normalize these components to —/j.b(N^ + N±), where iV-j = kp(l + I) 3 ^ 2 /(6n 2 ), and N± — 
fc3(l-/) 3 / 2 /(67r 2 ). 

The proximity effects can also be examined through the local DOS, N(y, e), given by, 

N( V , e) = - J2iKr(y) + <i(y)]f'(z - *0 + Kt (v) + tiMfte + (2-22) 

n 

1 



where /' = df/de. We will be concerned mainly with 
the DOS normalized to the DOS of a bulk (unpolarized) 
normal metal, Dn(0) — k 3 F / {2tt 2 E p) . 



III. RESULTS 

In this section we present our results, obtained self con- 
sistently as explained above and in previou o 18 ! 34 ' 35 work. 
We have assumed a coherence length fci?£o = 100. We 
will choose a geometry in which the layers are relatively 
thick: kpd s = D$ = 200 (that is, two coherence lengths) 
and Dp 1 = Dp 2 = 250. These values ensure that the 
sample will be overall superconducting at temperatures 
up to about 1/3 of the transition temperature T c , of a 
pure S bulk sample. This was not the case for the smaller 
values used in Ref. [l_8j where the condensation energy was 
quite small (see e.g. figure 6 in Ref. UK). This allows us 
to study the temperature dependence of the quantities 
involved over a broad range. The most important pa- 
rameters are the angle a and the magnet strength /. We 
will vary a in its full range between and tt and give re- 
sults for the values of / of 0.25, 0.5, and unity. No triplet 
amplitudes arise at I = (when magnetism is absent). 
In the results presented we have Zb = 0, when proximity 
effects arc in general maximized. 

In Fig. [2] we present comprehensive results for the real 
parts of fo(y,t) and f±(y,t), which we denote simply as 



fo{y,t) and fi(y,t) respectively. These are plotted in 
terms of the dimensionless variable Y = kpy. The am- 
plitudes are normalized to the value of the usual singlet 
amplitude in a pure bulk S sample. The temperature 
is set to zero in this figure. In the main plots, half of 
the S region and a portion (three fifths) of the left (F±) 
region are included. The corresponding portion on the 
F2 side can be inferred from the geometry and symme- 
try considerations. Results are plotted at three values 
of / and at a number of finite times r = u)p>t between 
0.4 and 8 as indicated in the legends. We have verified 
that at t = the computed triplet amplitudes vanish 
identically, in agreement with the Pauli principle. This 
is true, however, only when the calculation is performed 
to self-consistency: non-self consistent results invariably 
violate the Pauli principle near the interface. The results 
for /o are given at an angle a — while those for /1 are 
at a = 7r/2. At a = 0, f\ vanishes identically since the 
z component of the total spin is then a good quantum 
number. At a — n/2, and short time scales, the spatial 
dependences of the two triplet components coincide, al- 
beit with different signs in the two magnet regions, due 
to the magnetization vectors having equal projections on 
the x and z axes. At longer times, when the triplet am- 
plitudes extend throughout the S layer and couple the 
two magnets, /o and f\ deviate from one another. The 
insets in each panel amplify and clarify the region near 
the interfaces. 
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FIG. 2: The real parts, fo and /i, of the triplet pair ampli- 
tudes fo and fi (Eqns. (|2, 18(1 1. plotted as a function of position 
(in terms of Y = kpy) for three values of I at different times 
t = unt indicated in the legends of the top panels. These 
quantities are normalized to the value of the singlet pairing 
amplitude in a bulk S material. The main plots show half of 
the S region (right side of the vertical dashed line), and part 
of the Fi region. The insets are blow-ups of the region near 
the interface. The angle a is zero in the left panel and 7r/2 in 
the right panel. 
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FIG. 4: Maximum absolute values of fo and fi (see text) as a 
function of dimensionless time r at I — 1. In the top panel we 
consider fo at both a = and a = tv/2 while in the bottom 
panel we consider fi at a = tv/2. 
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FIG. 3: Penetration depths for the triplet amplitudes (see 
Eq. (|3.ip ). plotted as a function of r, as calculated from /o and 
fi in both the S and F regions, for the values of I indicated 
and the same angles as in Fig. [2] 



On the F side, both amplitudes peak very near the 
interface and then decay in an oscillatory manner, rem- 
iniscent of the behavior of the usual pair amplitude. 
Although the height of the first peak does not depend 
strongly on /, the subsequent decay in the F material is 
faster for larger values of /. This can be attributed to a 
decreased overall proximity effect: here we have assumed 
that at I = there would be no mismatch between the 
Fermi surface wavevectors of the two materials, implying 
that as / increases the mismatch between either the up 
or the down Fermi wavevectors k-\ and fcj, on the F side, 
and that in the S side increases. The location of this first 
peak depends very clearly on /, its distance to the inter- 
face decreasing as approximately 1 // consistent with the 
general rule that the oscillatory spatial dependences on 
the F side are determined by the inverse of fe-f — ki . The 
height of the first peak depends strongly on time and is 
maximum at times r of about 2tt. It is quite obvious that 
at intermediate values of I the penetration of the triplet 
correlation into the F material is rather long ranged. 

On the superconducting side the behavior is quite dif- 
ferent: the triplet correlations penetrate into the S mate- 
rial over a distance that rather quickly reaches two corre- 
lation lengths and then of course saturates at the sample 
size, without signs of decaying in time at these length 
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FIG. 5: Imaginary parts, fo(y,t) and fi{y,t) of the complex 
triplet amplitudes fo(y,t) and fi(y,t). These quantities are 
normalized and plotted exactly as their corresponding real 
parts are in Fig. [2j except that here we consider only the 
1 = 1 case, and both plots are for a = tt/2. As in Fig. [2] 
the main plot shows the behavior over an extended region. 
There are now two insets to each main plot, each showing the 
detailed behavior near the interface itself, on either the S or 
F side. 



scales. Furthermore this effect now increases sharply 
with / and is maximal in the half metallic case. Thus, 
the magnets act as sources, so to speak, of triplet correla- 
tions that enter the S material and this effect is stronger 
when / is larger. 

It is instructive to extract characteristic penetration 
lengths ii from the above data using the definition, 



fdy\My,t)\ 
max|/i(y,i) 



i = 0,l, 



(3.1) 



where the integration is either over the S or the F re- 
gion. In Fig. [3j the top two panels show the penetration 
lengths for the F material, at three values of / and the 
same values of a as for Fig. [2l The results are very sim- 
ilar whether they are calculated from the results for /o 
or from those for f\. The penetration length at constant 
time decreases with / as already noted and shows signs 
of saturating with time at a value which for / = 0.25 ap- 
proaches that of the superconducting coherence length. 



On the S side (bottom panels) the situation is very dif- 
ferent: the results for /o and /i are now clearly dissimilar 
with the penetration length for the former quantity being 
(for cases shown here) the larger one. This arises from the 
geometry and magnetization projections of each F layer 
on the x axis, which are in opposite directions, forcing 
the triplet fi to possess a node at the center of the tri- 
layer. No such requirement exists for f$, as it is spatially 
symmetric. Except for the case of fi at I = 1, we see 
no sign of saturation. In fact, the maximum value of r 
displayed here corresponds to the case in which the entire 
intrinsically singlet superconductor layer, two coherence 
lengths thick and sandwiched between two magnets, is 
wholly pervaded by induced triplet correlations. 

It is also of interest to consider the variation of the spa- 
tial maximum values of /o(y, t) and /i(y, t) with time. In 
Fig. 0] we show, for each time, the largest value of these 
quantities, which typically is attained near the interface, 
in either the F or S regions. By "maximum" value we 
mean the maximum of |/o| and not to be confused 
with the absolute value of the complex quantities |fo| 
or | f 1 1 . In this figure, the magnets are half metallic, 
1=1. In the top panel, we plot the results for /o at 
both a = and a = ir/2. We see that at earlier times, 
the value fo(y,t) at its peak just inside the F region (see 
Fig. [2 bottom left panel) exceeds the maximum value of 
this quantity in S. At longer times, however, there is a 
crossover as the size of the peaks in F decreases rather 
sharply, as explained above, while the size of the ampli- 
tude in S decreases only slowly, as the triplet correlations 
fill the S layer. It is apparent from careful examination of 
the / < 1 panels in Fig.[2J that this crossover does not oc- 
cur for smaller values of I except possibly on a time scale 
much longer than that considered here. In the bottom 
panel, we present a similar study of /i, this time of course 
only at a = tt/2 since this quantity vanishes identically 
for collinear magnetizations. The results are clearly very 
similar except that the results in the F region appear 
to saturate and do not decrease at long times. This is 
consistent with the earlier discussion where we saw that 
|/o| and |/i| overlap at a — 7r/2, except at sufficiently 
long times. In all cases the maximum value of the quan- 
tity plotted crests near r = 6 in agreement with previous 
remarks. 

All of the above results have been given in terms of the 
real parts, /o and /i of the complex amplitudes fo and 
fi . The behavior of the corresponding imaginary parts is 
qualitatively very similar and thus we will present only 
two examples, in Fig. [5] We denote these imaginary parts 
by fo(y,t) and fi(y,t) respectively. In Fig.[S]we consider 
only the case 1=1 (compare with Fig. [2]) and a = tt/2. 
As in Fig. [3 the main plots include an extended region 
near one of the interfaces and the insets are close views 
of the interface itself, in this case one of the insets shows 
a more detailed view of the S side. On the F side, the 
behavior is reminiscent to that of the real parts, except 
that the very prominent peak seen in the real parts right 
at the interface is absent for the imaginary parts. On the 
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FIG. 6: Magnetic moment component m x (see Eq. (|2.21|l '). 
normalized as explained in the text, plotted vs position (at 
1 = 1) for several values of a. The main plot shows the 
behavior near the interface (vertical dashed line), while the 
inset covers the whole sample. 



S side, the sign is now initially negative and it changes to 
positive at t of order unity. No such change was observed 
for the real parts. At longer times, the imaginary part of 
the triplet correlations also eventually penetrates several 
correlation lengths into the S sample, just as the real part 
does. 

In the next two figures, we explore the reverse prox- 
imity effect (the spreading of the magnetism into the S 
layer) as a function of a. This is best done by consider- 
ing separately the two components of the local magnetic 
moment vector. First, in Fig. [6] we consider the x com- 
ponent m x (see Eq. (|2.2ip ) normalized to the absolute 
value of its bulk value in a pure F material. The results 
in this figure are for half metallic magnets, (I — I). In 
the main plot of the figure we display the value of m x in 
the region very near an interface for several values of a. 
Of course, m x vanishes at a = 0. At other values of a 
it is large in the F material and it penetrates into S in 
an oscillatory way that is quite reminiscent of the corre- 
sponding penetration of the superconducting correlations 
into S. We see that the period of the spatial oscillations of 
the magnetization is independent of the angle a between 
the two F layer magnetizations. Another discernible fea- 
ture is that m x dampens out over relatively short length 
scales, consistent with past work*2i The inset shows the 
overall behavior of m x in the entire sample, demonstrat- 
ing also the opposite signs between in the two magnets 
in accordance with Fig. [TJ 

Similarly, in Fig. [7] we display the z-component of m z 
(see Eq. (|2.20p ). normalized in the same way as m x , and 
for the same values of a but including now three different 
values of /. Again, the main plots display the behavior 
near the interface while the insets are for the entire sam- 
ple. One can see here that the reverse proximity effect 
is very weak at small I and largest in the half metallic 



0.04 



-0.02 



N 0.04 



0.04 



0.00 




-100 -98 



-94 -92 -90 



Y 



FIG. 7: Normalized (see text) magnetic moment component 
m z (Eq. Q2,20p ) plotted vs. dimensionless position for three 
values of I. Again, the main plot is the behavior near the 
interface, while the insets cover the whole sample. It is evident 
that m z points in the same direction for both magnets 



case. The magnetic moment oscillates in the supercon- 
ductor with a period that is independent of the direc- 
tion and magnitude of the mutual magnetization in the 
F layers. The observed trends in m z hold also for the m x 
component. 

We next study the energy dependence of the single 
particle quasiparticle spectrum by considering the lo- 
cal density of states (local DOS), N(y,e), as defined in 
Eqn. (j2~22"l) . In Fig. we consider the local DOS, in- 
tegrated either over either the entire S or the entire F 
region, and normalized to its bulk value on a sample of 
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FIG. 8: Normalized local DOS from Eq. (|2.22|l integrated over 
(see text) over either the F region (left panels) or the S region 
(right panels) for three values of / and several relative mag- 
netization orientations. The temperature is at T = 0.01T C . 



the S material in its normal (non-superconducting) state. 
The results are displayed for three values of / and sev- 
eral values of a. The results reflect and confirm what we 
already have found out from analyzing the triplet ampli- 
tudes. On the F side, the proximity effect increases the 
correlations somewhat with a, as a reduction in quasipar- 
ticle states emerges for low energies due to increased cor- 
relations as the relative magnetizations become increas- 
ingly antiparallel. As / increases, the proximity effects 
weaken. Indeed, at I — 1 the DOS is nearly flat except 
at the larger values of a, suggesting that in the absence 
of a down spin band there may be a contribution from a 
possible triplet presence. On the superconducting side, 
the situation is somewhat different: the results never ap- 
proach a limit where the local DOS would look similar 
to that of a bulk superconductor. Even at / = 1, there is 
never a gap and even when that situation is approached 
at larger a, when the magnetizations in the two F lay- 
ers are antiparallel, the shape of the DOS curve does not 
resemble the signature bulk superconductor result. The 
profound difference between the parallel (a = 0) and an- 
tiparallel (a = 7r) cases is consistent in every respect with 
that previously^ found. 

The above results were all obtained in the low tem- 
perature limit. In the remaining figures we consider the 



FIG. 9: Temperature and position dependence of /o and /i at 
/ = 0.5 and r = 4. In the top panel a = and in the bottom 
panel a = n/2. The values of T/T c , where T c is the transition 
temperature of bulk S, are indicated. The main plots show 
a rather wide region near the interface (vertical dashed line) 
and the insets focus on the region very near the interface. 
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FIG. 10: The maximum values of the real parts of the triplet 
amplitudes at / = 0.5 and r = 4, as a function of temperature 
and the same values of a as in Fig. [9] The inset shows the 
corresponding peak values of the ordinary singlet amplitude, 
h = A(»)/ S (Eq. H23}). 
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temperature dependence. In Fig. [9] we plot directly the 
spatial behavior of /o and /i over a broad range of tem- 
peratures. As in previous figures for the triplet ampli- 
tudes, the main plot shows the behavior over a relatively 
extended region of the sample and the insets magnify the 
region near interface. In both cases we have taken / = 0.5 
and r = 4 while a — for /o and a = 7r/2 for fx. It can 
be inferred that at low T the temperature dependence 
of the triplet amplitudes is weak, while as T increases, 
the closer spacing in temperatures shown suggest that the 
correlations become destroyed at a much more rapid rate. 
In comparing the two panels, it is seen that for much of 
the temperature range, |/o| > |/i|, but at higher tem- 
peratures (T > 0.3T C ), f\ becomes the larger of the two. 
The temperature dependence of the characteristic pene- 
tration depths is weak. These observations are further 
exemplified in Fig. [101 where the peak values of |/o| and 
|/i| (also at r = 4) are shown as a function of T. These 
quantities are determined by calculating max{/ (y,T)} 
and max{/i(y, T)} throughout the structure for a given 
temperature. The inset depicts the corresponding peak 
values of the ordinary self- consistent equal-time singlet 
amplitude f 3 (y) = A(y)/g, (see Eq. (j!H)|) l For both the 
triplet and singlet behavior, there is a strong dependence 
on T as the temperature approaches the system's transi- 
tion temperature, which is about 0.32T C for our system. 
Technically, the determination of the self consistent am- 
plitudes is more difficult at higher T, when the number 
of iterations is in principle much higher. This increase in 
computational time can be reduced by up to an order of 
magnitude by taking as the initial spatial pair potential 
at a given T the result for the previously obtained next 
lower temperature times a T dependent factor derived 
from the linearized Ginzburg-Landau theory^. 

In accordance with what we have just seen, it is a natu- 
ral extension to study how thermal effects might destroy 
the spatial characteristics of singlet correlations through- 
out the structure. In Fig. [TTJ we therefore display /a as 
a function of Y, for several temperatures, and in which 
the relative magnetizations are collinear (a = 0) and at 
right angles (a = 7r/2). Remarkably, one can see that 
the temperature dependence of the triplet components is 
somewhat weaker that that of the standard singlet ampli- 
tude. One can also clearly see that, as indicated above, 
the penetration of the triplet amplitudes into the F ma- 
terial over a length scale that is clearly much longer than 
that of the singlet amplitude. This is again a strong in- 
dication that experimental tunneling results indicating 
long ranged penetration effects in F/S structures are in- 
deed evidence for induced triplet correlations. 



IV. CONCLUSIONS 

In this paper, we have presented a detailed study of 
induced time dependent (odd in time) triplet pairing 
correlations in clean planar FSF junctions consisting of 
an ordinary s-wave superconductor sandwiched between 




-200 -150 -100 -50 



Y 

FIG. 11: Temperature and position dependence of the or- 
dinary singlet pair amplitude (or A (y) /<?)), normalized to its 
value in bulk S material, for moderate magnetic strength, I = 
0.5. Top panel: the magnetizations are parallel (a = 0). Bot- 
tom panel: the magnetizations are perpendicular (a = n/2). 



two relatively thick ferromagnets whose magnetizations 
are misoriented with respect to each other by an angle 
a. Our microscopic formalism allowed us to investigate 
cases involving strong magnets, as well as atomic scale 
phenomenon, two things not possible in the widely used 
quasiclassical approaches. We have obtained results as a 
function of a, time, the strength I of the ferromagnets, 
and the temperature. We have presented results for the 
spatial behavior of the time-dependent triplet pair ampli- 
tudes, /o and fx, and for the corresponding penetration 
lengths extracted from them. We have found that these 
triplet correlations arc indeed induced via the proxim- 
ity effect, that they completely pervade even a super- 
conductor several coherence lengths thick, and also sub- 
stantially penetrate the ferromagnetic layers. These re- 
sults have clear implications for the experimental work in 
which long range proximity effects in SF nanostructures 
have been reported, effects that have been speculated to 
be due to the existence of some kind of triplet pairing. 
Our calculations, in which the time dependence is stud- 
ied from the Heisenberg picture, emphasize the need for 
full self-consistency of the solutions, without which we 
hnd that the Pauli principle is violated. 

We have also considered the reverse proximity effect, 
which is of particular interest in this case due to the pres- 
ence of two components of the magnetization, and we also 
have given results for the experimentally measurable lo- 



12 



cal density of states, which revealed clear subgap energy 
signatures as a function of a and /. We have studied 
the temperature dependence of the triplet amplitudes (as 
well as the ordinary singlet amplitude) and found that 
the temperature dependence of the penetration depths 
associated with these triplet amplitudes is weak: these 
lengths remain large all the way up to the vicinity of the 
transition temperature of the system. This bodes well 
for further experimental observations and verification in 
these clean systems. 
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